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Abstract: SOFTSUSY is a program which accurately calculates the spectrum of superpar- 
ticles in the CP-conserving Minimal Supersymmetric Standard Model (MSSM), with a full 
flavour mixing structure. The program solves the renormalisation group equations with 
theoretical constraints on soft supersymmetry breaking terms provided by the user. Weak- 
scale gauge coupling and fermion mass data (including one-loop finite MSSM corrections) 
are used as a boundary condition, as well as successful radiative electroweak symmetry 
breaking. The program can also calculate a measure of fine-tuning. The program structure 
has been designed to easily generalise to extensions of the MSSM. This article serves as a 
self-contained guide to prospective users, and indicates the conventions and approximations 
used. 



Keywords: ^particle, MSSM . 



Contents 



i 



i 



Introduction 

The Nature of the Physical Problem 



1.2 The Program 



1.3 Aims and Layout 



MSSM Parameters 
|2.l| Super symmetric Parameters 
SUSY Breaking Parameters 
Tree-Level Masses 



2.2 



2.3 



3.2 



Calculation 
|3l| Below Mz 

Initial Estimate 
Gauge and Yukawa Couplings 
MSSM Renormalisation 
Electroweak Symmetry Breaking 



3.3 



3.4 



3.5 



3.6 MSSM Spectrum 



3.7| Fine Tuning 



1^ Running SOFTSUSY 

lT| Input file 



Sample Output 



B.l Problem flags 



[^ . Sample Program 

|D| . Switches and Constants 



Object Structure 

Linear Algebra 
General Structure 
RGE Class 



E.l 



E.2 



E.3 



EL4| QedQcd Class 

E.5 MssmSusy Class 

R6| Sof tParsMssm Class 

E.7 MssmSof tsusy Class 

E.8 MssmSof tsusyAltEwsb Class 

E.9 FlavourMssmSof tsusy Class 







m 

10 

n 

12 

13 
11 

15 
P 

19 
22 
22 

n 

22 

m 



- 1 - 



1. Introduction 



The Minimal Supersymmetric Standard Model (MSSM) provides an attractive weak-scale 
extension to the Standard Model. As well as solving the gauge hierarchy problem, it can be 
motivated by more fundamental models such as various string theories or supersymmetric 
grand unified theories. The MSSM provides a rich and complicated phenomenology. It 
predicts many states extra to the Standard Model (sparticles) and their indirect empirical 
effects and direct detection are vital for verification of the MSSM. Models that are more 
fundamental than the MSSM can provide stringent constraints upon the way supersymme- 
try (SUSY) is broken, with important implications for the spectrum which in turn affects 
the signatures available in experiments. It is therefore desirable to construct a calculational 
tool which may provide a spectrum and couplings of the MSSM sparticles so that studies 
of the capabilities of colliders, extraction of high scale parameters (if a signal is observed) 
and studies of constraints on the models are enabled. We present such a tool (SOFTSUSY) 
in this article. 

1.1 The Nature of the Physical Problem 

The determination of sparticle masses and couplings of SUSY particles in the R-parity con- 
serving MSSM is the basic problem. Low energy data on Standard Model fermion masses, 
gauge couplings and electroweak boson masses are to be used as a constraint. SUSY radia- 
tive corrections to these inputs from sparticle loops depend upon the sparticle spectrum, 
and must be calculated. Theoretical constraints on the SUSY breaking parameters from 
an underlying theory are often imposed at a high renormalisation scale, perhaps resulting 
from a super gravity or string theory. Often, the theoretical constraints drastically reduce 
the number of free parameters in the SUSY breaking sector (which numbers over 100 in 
the unconstrained case). These constraints then make phenomenological analysis tractable 
by reducing the dimensionality of parameter space sufficiently so that parameter scans 
over a significant volume of parameter space are possible. Finally, the MSSM parameters 
must also be consistent with a minimum in the Higgs potential which leads to the observed 
electroweak boson masses. 

This problem has been addressed many times before in the literature (see for exam- 
ple m, ^, ^, ^ P), with varying degrees of accuracy in each part of the calculation. It is our 
purpose here to provide a tool which will solve the problem with a high accuracy, including 
state-of-the-art corrections. Similar problems in the context of MSSM extensions^ have 
also been studied. In anticipation of new forms of SUSY breaking constraints and new 
MSSM extensions, we designed the tool to be flexible and easily extended. 



1.2 The Program 

SOFTSUSY has been written in object-oriented C++ but users may use an executable pro- 
gram with input either in the SUSY Les Houches Accord [29| format or from command- line 
arguments. For users wishing to call SOFTSUSY from their own programs, the user interface 



^By MSSM extension, we mean an extension applicable near the weak scale. 
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is designed to be C-like to aid users that are unfamiliar with object orientation. Accuracy 
and ease of generalisation have taken priority over running speed in the design. For exam- 
ple, full three family mass and Yukawa matrices may be employed, rather than the more 
usual dominant third family approximation. The publicly released codes ISASUGRA (which 
comprises part of the ISA JET package [^) and SUSPECT [Q] use the dominant third family 
approximation, for example. The full three-family choice slows the renormalisation group 
evolution significantly, but will facilitate studies of sparticle or quark mixing. The running 
time is not foreseen as a bottleneck because it is a matter of a couple of seconds on a 
modern PC, and will certainly be negligible compared to any Monte-Carlo simulation of 
sparticle production and decay in colliders. It is possible for the user to specify their own 
high scale boundary conditions for the soft SUSY breaking parameters without having to 
change the SOFTSUSY code. For the convenience of most users however, the most commonly 
used high-scale boundary conditions are included in the package. 

The code can be freely obtained from the SOFTSUSY web-page, which currently resides 
at the address 



http:/ /projects. hepforge.org/softsusy/ 

Installation instructions and more detailed technical documentation of the code may also 
be found there. 

SOFTSUSY is a tool whose output could be used for studies of MSSM sparticle searches Q 
by using event generators such as HERWIG 0, or other more theoretical or astrophysical 
studies. For a review of SUSY tools on offer (which may use the output from SOFTSUSY), 
see Ref. [|]. 

1.3 Aims and Layout 

The main aims of this article are to provide a manual for the use of SOFTSUSY, to describe 
the approximations employed and to detail the notation used in order to allow for user 
generalisation. There have been other articles published on the comparison of the calcu- 
lation in SOFTSUSY with those of other codes and so we decline from including such 
information here. 

The rest of this paper proceeds as follows: the relevant MSSM parameters are presented 
in sec. |2|. The approximations employed are noted in sec. ^, but brevity requires that they 
are not explicit. However, a reference is given so that the precise formulae utilised may 
be obtained in each case. The algorithm of the calculation is also outlined. Technical 
information related to running and extending the program is placed in appendices. A 
description of how to run the command-line interface is given in appendix |^, including 
information on the input-file. Sample output from one such run is displayed and explained 
in appendix Appendix |^ gives a more technical example of a main program, useful if 
the user wants to call SOFTSUSY from his or her main program. The use of switches and 
constants is explained in appendix Finally, in appendix ^ a description of the relevant 
objects and their relation to each other is presented. 
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2. MSSM Parameters 



In this section, we introduce the MSSM parameters in the SOFTSUSY conventions. Trans- 
lations to the actual variable names used in the source code are shown in appendix ^ 

2.1 Supersymmetric Parameters 

The chiral superfields of the MSSM have the following Gsm = SU{2,)c x SU{2)l x U{1)y 
quantum numbers 

L:(l,2,-i), ^: (1,1,1), Q: (3,2,1), [7: (3,1,-^), 

5:(3,l,i), i^i : (1,2,-i), H2:{1,2,^). (2.1) 

Then, the superpotential is written as 



{YE)ijL'iHlEj + {Yr>)ijQf^H1D^,, + {Yu)ijQrmUj^ + ^iH^Hl (2.2 



Throughout this section, we denote an SU{3) colour index of the fundamental represen- 
tation by x,y,z = 1,2,3. The SU{2)i fundamental representation indices are denoted by 
a, 6, c = 1, 2 and the generation indices hy i, j, k = 1,2,3. e^b = e"'' is the totally antisym- 
metric tensor, with ei2 = 1. Note that the sign of fi is identical to the one in ISASUGRA 
but is in the opposite convention to ref. @. Presently, real Yukawa couplings only are 
included. All MSSM running parameters are in the DR scheme. The Higgs vacuum expec- 
tation values (VEVs) are {Hf) = Vi/\/2 and tan/3 = V2/V1. gi are the MSSM DR gauge 
couplings and gi is defined in the Grand Unified normalisation gi = y^5/3g' , where g' is 
the Standard Model hypercharge gauge coupling. Elements of fermion mass matrices are 
given by 

{mu)ij = -j={Yu)ijV2, imd,e)ij = -^{YD,E)ijVi (2.3) 
for the up quark, down quark and charged lepton matrices respectively. 
2.2 SUSY Breaking Parameters 

We now tabulate the notation of the soft SUSY breaking parameters. The trilinear scalar 
interaction potential is 



V3 = eab 



QT,{UA)ijUj,^Hl + Qf^iDA^dj^^H^ + Ll {EA)ifej^m + H.c. , (2.4) 



where fields with a tilde are the scalar components of the superfield with the identical 
capital letter. Also, 

{Au,D,E)ij = {Ua, Da, EA)ij / {Yu,D,E)ij (2.5) 

(no summation on are often referred to in the literature as soft j4-parameters. 
The scalar bilinear SUSY breaking terms are contained in the potential 



)ijUj^ -\- (ij ijn^ijdj^ -\- ei{mg)ijej 



u\{mi)iju*.^ + d1{m%^d*^ + ei(m|)ije* + eah{vallilR\ + R.c). (2.6) 
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For a comparison of these conventions with other popular ones in the hterature, see Table |. 
In the table, we compare the SOFTSUSY conventions with the SUSY Les Houches Accord p9| 
and ref. (Martin and Vaughn). 

Writingthebinoas6,u;^=i'2,3as the sqptsusy SLHA Martin and Vaughn 
unbroken-SU(2)z, gaugmos and yC/,D,s '(j1J:d:ey 

as the gluinos, the gaugino mass terms j, 

^ ■ . ■ T ■ Ua,-L>A,^A J-U,D,E "'U,D,E 

are contained m the Lagrangian r, r, r, 

CG = -[Mibb + M2w''w^ + M3~g''~g'')+h.c. "^J,- "^h/ <J,e 

(2.7) fJ' 

2.3 Tree-Level Masses 

Table 1: Comparison of conventions between 
Here we suppress any gauge indices and SOFTSUSY and the literature, 
follow the notation of ref. closely. The 

Lagrangian contains the neutralino mass matrix as —^ip^'^Ad^oW + h.c, where ip^ = {—ib, 
—iw^, hi, ^2)^ and 



/Ml -MzCf3Sw Mzspsw \ 

M2 Mzcpcw —Mzsfjcw 

-MzCfssw MzCfscw -^u 
\ Mzspsw -Mzspcw -M / 



(2. 



We use s and c for sine and cosine, so that sp = sin/3, = cos j3 and sw{cw) is the sine 
(cosine) of the weak mixing angle. The 4 by 4 neutralino mixing matrix is an orthogonal 
matrix O with real entries, such that M^oO is diagonal. The neutralinos Xi* are defined 
such that their absolute masses increase with increasing i. Some of their mass values can 
be negative. 

We make the identification w'^ = {w^ =F iw^) /V^ for the charged winos and , /i^ for 
the charged higgsinos. The Lagrangian contains the chargino mass matrix as + 
h.c, where = {—iw~^, h'^)'^ 1 ^~ ~ h^)^ and 



M2 V2MwsA ,^ 
V2MwCf3 ^^ J' 



This matrix is then diagonalised by 2 dimensional rotations through angles 0^,9 ft in the 
following manner: 



mt.. 



where mj, could be negative, with the mass parameter of the lightest chargino being in 
the top left hand corner. 

At tree level the gluino mass, m^, is given by M3. 

Strong upper bounds upon the intergenerational scalar mixing exist ||To[| and in the 
following we assume that such mixings are negligible. The tree-level squark and slepton 
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mass squared values for the family i are found by diagonalising the following mass matrices 
Ad'j defined in the {fiL-,fiR}^ basis: 

(i^%)ii + ^li + i\-hwWhw mui{{Au)ii- n cot 13) \ 
\ ruui {{Au)ii - /xcot /3) (jnDu + ml. + |s^M|c2^ j 



+ m| - (i - is^)M|c2/3 mdi{{AD)ii- lJitSMl3) \ 
{{AD)ii - /xtan/3) (mj)ji + m^. - is^M|c2/3 / 

/ (m|)ii + m^^ - (i - s^)M|c2/3 me, ((^B)n - ^tan/3) ^ 
I met {(AEh - Aitan/3) (m?)n + m^, - s^M|c2/3 - 



(2.11) 
(2.12) 



(2.13) 



mf,ef are the mass and electric charge of fermion / respectively. The mixing of the first 
two families is suppressed by a small fermion mass, which we approximate to zero. The 
sfermion mass eigenstates are given by 



mj^ 
mj^ 




(2.14) 



where is the cosine of the sfermion mixing angle, cos Of, and Sf the sine. 9f are set in the 
convention that the two mass eigenstates are in no particular order and 6f e [— 7r/4, 7r/4]. 
The sneutrinos of one family are not mixed and their masses are given by 

m% = {m\)u + ^M|c2^. (2.15) 

The CP-even gauge eigenstates (i?°, i?^) are rotated by the angle a into the mass eigen- 
states [H^ h^) as follows, 

m^o < mjjo by definition, and a G [— 7r/4, 37r/4]. The CP-odd and charged Higgs masses 
are 

m\o = m|(tan (3 + cot /?), mjj± = m\o + (2-17) 




at tree level. 



3. Calculation 

We now show the algorithm used to perform the calculation. Standard Model parameters 
(fermion and gauge boson masses, the fine structure constant a{Mz)-, the Fermi constant 
from muon decay G^, and a2,{Mz)) are used as constraints. The soft SUSY breaking 
parameters and the supcrpotential parameter /x are then the free parameters. However, in 
what follows, is constrained by Mz and tan/3 is traded for m^ as an input parameter. 
Therefore, the total list of unconstrained input parameters is: any fundamental soft SUSY 
breaking breaking parameters (except m|), tan /3 and the sign of /i. First we describe the 
evolution of the low-energy Standard Model input parameters below Mz, then detail the 
rest of the algorithm. 
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3.1 Below Mz 



a{Mz), as{Mz) are first evolved to 1 GeV using 3 loop QCD and 1 loop QED in|, y, |13| 
with step-function decoupling of fermions at their running masses. We have checked that 
the contribution from 2-loop matching |14| is negligible; the effect of 3- loop terms in the 



renormalisation group equations is an order of magnitude larger. Then, the two gauge 
couplings and all Standard Model fermion masses except the top mass are run to Mz- The 
f3 functions of fermion masses are taken to be zero at renormalisation scales below their 
running masses. The parameters at Mz are used as the low energy boundary condition in 
the rest of the evolution. 

3.2 Initial Estimate 

The algorithm proceeds via the iterative method, and therefore an approximate initial 
guess of MSSM parameters is required. For this, the third family DR Yukawa couplings 
are approximated by 

h.(Q) = hAQ) = ^^^i^. (3.1) 

V sm p V cos p 

where v = 246.22 GeV is the Standard Model Higgs VEV and Q = mt{mt) is the 
renormalisation scale. The AdS values of fermion masses are used for this initial esti- 
mate. The fermion masses and at the top mass are obtained by evolving the pre- 
viously obtained fermion masses and gauge couplings from Mz to rrit (with the same 
accuracy). The electroweak gauge couplings are estimated by ai(Mz) = 5a(M^)/(3c^), 
a2{Mz) = a{Mz)/s^i. Here, sw is taken to be the on-shell value. These two gauge cou- 
plings are then evolved to mt with 1-loop Standard Model (3 functions, including the effect 
of a light higgs (without decoupling it). In this initial guess, no SUSY threshold effects 
are calculated. The gauge and Yukawa couplings are then evolved to the unification scale 
Mx with the one-loop MSSM (3 functions, where the user-supplied boundary condition on 
the soft terms is applied. Also, fi{Mx) = sgn(/i) x 1 GeV and m3{Mx) = are imposed. 
These initial values are irrelevant; they are overwritten on the next iteration by more real- 
istic boundary conditions. /i(Mx) is set to be of the correct sign because its sign does not 
change through renormalisation. 

The whole system of MSSM soft parameters and SUSY couplings is then evolved 
to 1-loop order to Mz- At Mz, the tree- level electroweak symmetry breaking (EWSB) 
conditions are applied P to predict and 771,3. The masses and mixings of MSSM super 
particles are then calculated at tree-level order by using the SUSY parameters (and ms) 
calculated at Mz- The resulting set of MSSM parameters is then used as the initial guess 
for the iterative procedure described below. 

3.3 Gauge and Yukawa Couplings 

Figure m shows the iterative procedure, starting from the the top. The whole calculation is 
currently performed in the real full three family approximation, i.e. all Yukawa couplings 
are set to be real, but quark mixing is incorporated. First of all, the one-loop radiative 
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SUSY radiative corrections to gi{Mz), ht^b,T{Mz)- 



Run to Mx- Apply soft SUSY breaking boundary condition. 



convergence 



Run to Ms- 



EWSB, iterative solution of jj, 



Run to Mz- 



Calculate Higgs and sparticle pole masses. Run to Mz- 



Figure 1: Iterative algorithm used to calculate the SUSY spectrum. Each step (represented by a 
box) is detailed in the text. The initial step is the uppermost one. Ms is the scale at which the 
EWSB conditions are imposed, as discussed in the text. Mx is the scale at which the high energy 
SUSY breaking boundary conditions are imposed. 



corrections are applied to the gauge and third-family Yukawa couplings. For these, we 
rely heavily on ref. ^ by Bagger, Matchev, Pierce and Zhang (BMPZ)^. In the threshold 
corrections, we use running DR masses and parameters at the relevant scale, unless denoted 
otherwise in the text. mt{Mz) is calculated with 2-loop QCD |15] and the full one-loop 
supersymmetric contributions to mt{Mz), including logarithmic and finite contributions 
(eqs. (D.I6)-(D.I8) of BMPZ): 



MMz)'iPssM 



as{Mz) 



Sir 



(5 - 3L) + 



43 



ai{Mz){-0.538 + ^L-^L 



87r2' 



(3.2) 



where L = ln{m^{Mz)^/gg]^.j/M^). We denote the BPMZ corrections without the one-loop 
QCD part as Tif^^^^. These corrections are necessary because the region of valid EWSB 
is very sensitive to [0. To calculate mh{Mz)^^, we first calculate the Standard Model 



DR value from the MS one [15, 32 1 

mb{Mz)Wi = MMz)Wii^ 



a 



DR 



23a 



DR^ 



37r 



13gf 

727r2 ' I287r2 I1527r2 



+ 



(3.3) 



We then add the leading one-loop supersymmetric corrections 

mbiMz)WssM = mb{MzM/{l + A^^^^). 



(3.4) 



^Whenever a reference to an equation in BMPZ is made, it is understood that the sign of fi must be 
reversed. 
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The contributions to ^'susY included in full from eq. D.18 of BMPZ (neglecting the 
term proportional to e, since that is already included in the QED calculation of the SM 
mf,{Mz))- Both finite and leading logarithmic corrections are included. After the Standard 
Model MS bar of m,- is converted to the DR value via 

mAMzM = mAMz)Wi (l " J^ig! " ^l)) • (3-5) 

The full one-loop MSSM corrections from the appendix of BPMZ (aside from the photon 
contribution, since that has already been included in mr{Mz)g^j^i) are then used to correct 
m,(Mz)|i: _ _ 

mAMz&M = mr{Mz)^§{l + S,). (3.6) 



The one-loop DR values for mt{Mz), mb{Mz), mT-{Mz) are then substituted with the DR 



value of v{Q) into eq. ( |2.3D to calculate the third family DR Yukawa couplings at Mz- 
v{Q) is run to two-loops. The other diagonal elements of the Yukawa matrices are set by 



eq. ( p. 3D but with fermion masses replaced by the MS values. 

The default option is to perform the calculation in the dominant third-family approx- 
imation, where all elements of Yukawa matrices expect for the (3,3) elements are set to 
zero. There are also options described in appendix ^ for performing the calculation in the 
unmixed 3-family approximation or the fully-mixed 3-family case. If a flavour-mixing op- 
tion is chosen (see section [E.9D , the Yukawa couplings are then mixed using the "standard 
parameterisation" of the CKM matrix with CP-violating phase set either to zero or vr, 
whichever results in a positive entry for {Vckm)i3 (also known as Vub)'- 



VcKM 



^ C12C13 S12C13 psiz 

-S12C23 -PC12S23S13 C12C23 -PS12S23S13 S23C13 h (3-7) 

\ S12S23 -PC12C23S13 -C12S23 -P'Si2C23Si3 C23C13 



where Sij = sin^jj, Cij = cos Oij and p = ±1. Sign conventions are automatically chosen 
such that diagonal entries and entries above the diagonal in Vckm are positive. Note 
that Vckm is a member of 0(3), i.e. V^^j^ = V^j^j^j. While complex phase effects are 
obviously not taken into account in eq. |3.7| , it is hoped that the magnitudes of the main 
quark mixing effects will be. In fact, using the central values for 6ij given by the particle 



data group the magnitudes of all elements in Vckm in the first row and last column 
are exactly reproduced. Of the other entries, \Vcs\ is accurate at the 2xl0~^ level and 
\Vcd\ at the 4xl0~^ level: surely negligible for most practical purposes. \Vts\ is also quite 
accurate (to 1.2%), but it should be noted that \ Vtd\ is wrong by around 50%. Any flavour 
physics effects sensitive to \Vtd\ is therefore subject to this large uncertainty on its value 
issuing from SOFTSUSY^. The up (by default), or down Yukawa couplings at Mz are mixed 
in the weak eigenbasis via 

(yuY = V^km(Y'')Vckm, {Yd)' = Vckm{Y'')V^km (3.8) 



^It is hoped in the future to include complex phases in Yukawa matrices, sfermion soft mass squared 
terms and trilinear scalar couplings. 
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where the primed Yukawa matrix is in the weak eigenbasis and the unprimed is in the mass 
eigenbasis. 

Full one-loop corrections to gi{Mz) are included. The treatment of electroweak gauge 
couplings follows from appendix C of BMPZ, and includes: two- loop corrections from the 
top, electroweak boson and the lightest CP-even Higgs. The pole value of mt is used in 
the calculation of the W and Z self energy in order to calculate sin^^, since this is what 
is assumed in the two-loop corrections. We use the fine structure constant a{Mz)'^^^ , the 
Z-boson mass Mz and the Fermi decay constant as inputs. Mw is predicted from 
these inputs. Because the EWSB constraints tend to depend sensitively upon gi^2{Mz), 
accurate values for them are determined iteratively. An estimate of the DR value of Syy is 
used to yield a better estimate until the required accuracy is reached (usually within 3 or 
4 iterations). The QCD coupling is input as as{Mz)^^^ and is modified by gluino, squark 
and top loops as in eqs. (2), (3) of BMPZ in order to obtain the MSSM DR value. 

3.4 MSSM Renormalisation 

All soft breaking and SUSY parameters are then evolved to the scale 



Ms = ^mi^{Ms)mi^{Ms), (3. 



where [^] the scale dependence of the electroweak breaking conditions is smallest. Through- 
out the iteration described here, the renormalisation group evolution (RGE) employs three 
family, 2-loop MSSM (5 functions for the supersymmetric parameters tan/3 and the 
Higgs VEV parameter v are also run to two-loop order in the Feynman gauge, although 
the Higgs VEV RGE is missing terms C'(5'2,525i'9i)/(167r^)^ |2^. The program can 
be made much faster by switching certain two-loop terms off (scalar sfermion masses and 
trilinear couplings) as described in section j^. There is no step-function decoupling of spar- 
ticles: this is taken into account at leading logarithmic order in the radiative corrections 
previously calculated at Mz and in the calculation of the physical sparticle spectrum at 
Ms, described below. All (3 functions are real and include 3 family (and mixing) contribu- 
tions. If no flavour mixing is present in the model specified by the user, the 2-loop parts 
of the RGEs switch to the dominant 3rd-family version, where the lighter two families 
Yukawa couplings are neglected. 

3.5 Electroweak Symmetry Breaking 

The Higgs VEV parameter v{Ms) is set by: 

„.(M,) = 44a±^I%«lL, (3.10) 

where Mz is the pole Z mass and H^^ is the transverse Z self-energy. The full one-loop 
EWSB conditions at this scale are then employed to"^ calculate m3(Ms) and fi{Ms). fJ,{Ms) 

*Note that there is also an option to extract niHiiMsuSY) and niH^iMsu sy) from input fJ-{MsuSY) 
and myi(pole) values, see section [ ' 
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requires an iterative solution because the tadpoles depend upon the value of /i assumed. 
The symmetry breaking condition for /i can be phrased as 0] 



fJ-^ = tan 2/3 



tan /3 — cot /3 



M| , (3.11) 



where = m|^, — ti/vi, = + 5Rn^^(M^) is the running Z mass and ti are the 



tadpole contributions. The value of /x coming from the tree-level EWSB condition (eq. [3.11 



with 3f?n^^ = ti = 0) is utilised as an initial guess, then the one-loop contributions in the 
tadpoles and self-energy terms are added to provide a new value of fi{Ms). Two-loop 
terms 0{af), O^abOr), O(a^), 0{ahas), 0{atas), 0{a'^) and 0{ata},) are also included 



in the tadpoles |2l|, The tadpole corrections are then calculated using the new value 
of n{Ms) and the procedure is repeated until it converges to a given accuracy. m^lMs) is 
then determined by input the value of n{Ms) into the EWSB condition 

^l = ^(^k +ml+2f,'). (3.12) 



2 \ -Hi ' -H2 

The ensemble of MSSM parameters are then evolved using the (3 functions described above 
to the user supplied scale Mx- If gauge-unification has been specified as a boundary 
condition, the current estimate of Mx is revised to leading log order to provide a more 
accurate value upon the next iteration: 

M^^"" = Mxe«i(^^)-«2(^^) , (3.13) 

where primes denote derivatives calculated to 2-loop order. The user-supplied boundary 
conditions are then imposed upon the soft terms before the model is evolved back down 
to Ms- The super particle mass spectrum is determined at this scale. Because fx and ma 
are more scale independent at Ms as opposed to some other scale, the Higgs, neutralino 
and chargino masses also ought to be more scale independent by determining them at this 
scale. 



3.6 MSSM Spectrum 

In the following description of the approximations involved in the calculation of the super 
particle spectrum, it is implicit that where masses appear, their DR values are employed. 
The running value of swip) = e(//)/5f2(/^) is also employed. In loop corrections to spar- 
ticle masses, the Yukawa couplings of the first two families are set to zero, being highly 
suppressed compared to those of the third family. All sparticle masses are calculated with 
the full SUSY one-loop BPMZ corrections at the scale MsusY- Most sparticle masses are 
calculated at external momenta equal to their DR mass m{MsusY)- 

The physical gluino mass is calculated to full one-loop order as follows. The running 
parameters are evaluated at renormalisation scale ^ = MsusY and external momentum 
p = M^{ii) in the following corrections: 

= ifef f 15 + 9 In (^) - ^ B,{p, m„ ^)- 

\ q i=l 
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^ 1 ^ S2e, [Bo{p,mq,mq^,fj,) - Boip,mq,mq^, fi)] . (3.14) 

q=t,b ^^^^ J 

The Passarino-Veltman functions -Bo,i are given in appendix B of BMPZ. The physical 
gluino mass is then given by 

m-g = M^{Msusy) (1 + ^g{MsusY)) . (3.15) 

The gluino mass is allowed to be negative, as is the case in mAMSB, for example. Of 
course the kinematic mass is just the absolute value, and the phase may be rotated away, 
altering the phases of some of the Feynman rules. Negative masses for neutralinos can also 
be rotated away in this way. 

For mixed sparticles (sbottoms, staus, neutralinos and charginos), the external mo- 
mentum is equal to the DR value of the diagonalised mass. All mixing angles and matrices 
are defined such that they diagonalise the one-loop corrected mass matrix evaluated at the 
minimum DR mass eigenvalue. When flavour mixing is used, at tree-level, the full flavour 
structure is present in sfermion masses. Loop corrections are added only to the equal family 
pieces, not to intra-family mixing pieces as an approximation. All of the SOFTSUSY loop 
corrections also neglect intra-family mixing themselves. 

The pseudo-scalar Higgs mass m^o and CP-even Higgs masses m^o , rrifjo are deter- 
mined to full one-loop order as in eq. (E.6) of BMPZ in order to reduce their scale depen- 
dence, which can be large ||2^. The zero-momentum 0{al), 0{ahOir)t 0{a1), 0{ahas), 
0{atas), 0{a^), 0{atab) 2-loop corrections are also included in mjffi^m^o,mHO [|^, 



25| , Pq ], All one-loop corrections are included in the determination of the charged Higgs 



pole mass. Every Higgs mass is determined at an external momentum scale equal to its 
DR mass. 

Finally, the running MSSM parameters are evolved back down to Mz- The whole 
process is iterated as shown in flgure ||, until the DR sparticle masses evaluated at Ms all 
converge to better than the desired fractional accuracy (TOLERANCE), which may be set by 
the user in the main program or input flle. 

3.7 Fine Tuning 

We now detail the fine-tuning calculation. As lower bounds on super partner masses are 
pushed up by colliders, and m/^j niay be forced to be much larger than Mz if they are 
related to the other super particle masses, as is the case for example in the case of minimal 
super gravity. If we re-phrase eq. ( 3.11| ) as 



M| = -2^i^ + tan 2j3 m^^ tan /? - m^^ cot j3 , (3.16) 

we see that the terms on the right-hand side must have some degree of cancellation in order 
to reproduce the observed value of Mz- But ^ has a different origin to the SUSY breaking 
parameters and the balancing appears unnatural. Various measures have been proposed 



in order to quantify the apparent cancellation, for example ref.s |27, The definition of 
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naturalness Ca of a 'fundamental' parameter a employed here is [25 



a In M| 



5 In a 



(3.17) 



From a choice of a set of fundamental parameters defined at the scale Mx'- {fli}; the 
fine-tuning of a particular model is defined to be c = max(C(j). {oj} are any parameters 
in the user supplied boundary condition on the soft supersymmetry breaking parameters 
augmented by /if (Mx), /x(Mx) and m^{Mx)- The derivatives in eq. ( |3.17] ) are calculated 
by numerically finding the derivative of M^°'*^ = Mz + 3^n|^^(M^) in eq. ( |3.11 ). The 



input parameters are changed slightly (one by one), then the MSSM parameter ensemble 
is run from Mx to Mg where the sparticle mass spectrum is determined along with the 
corresponding MS Higgs VEV parameter = vf+V2- First of all, tan (3{Ms) is determined 
by inverting eq. ( 3.12| ) and the resulting value is utilised in a version of eq. ( ^.11 ) inverted to 



give M^^^ in terms of the other parameters. The resulting value of M^"'^ is the prediction 
for the new changed input parameters, and its derivative is determined by examining its 
behaviour as the initial changes in input parameters tend to zero. 

A. Running SOFTSUSY 

A main program is supplied which produces an executable called sof tpoint . x, is included in 
the SOFTSUSY distribution. For the calculation of the spectrum of single points in parameter 
space, we recommend the SLHA input/output option. The user must modify a file 
(e.g. lesHouchesInput, as provided in the standard distribution) that specifies the input 
parameters. The user may then run the code with 

. /sof tpoint . X leshouclies < lesHouchesInput 

In this case, the output will also be in SLHA format. Such output can be used for input 



into other programs which subscribe to the accord, such as SDECAY [34|, PYTHIA ||35[] (for 
simulating sparticle production and decays at colliders) or mirOMEGAs ||3^ (for calculating 
the relic density of neutralinos, 5 — > 57 and cy), for example. For further details on 
the necessary format of the input file, see ref. [^] and appendix [A.lJ . 

If the user desires to quickly run a single parameter point in AMSB, mSUGRA or 
GMSB parameter space, but does not wish to use the SLHA, the following options are 
available: 

. /sof tpoint . X sugra <mO> <ml2> <aO> <tanb> <mgut> <sgnMu> 

. /sof tpoint . X amsb <mO> <m32> <tanb> <mgut> <sgnMu> 

. /sof tpoint . X gmsb <n5> <mMess> <lajiibda> <tcLnb> <sgnMu> 

Bracketed entries should be replaced by the desired numerical values, (in GeV if they are 
dimensionful) . The program will provide output from one point in mSUGRA, mAMSB or 
GMSB. If <mgut> is specified as unified on input, SOFTSUSY will determine mgut to be the 
scale that Qi^Mgut) = g2iMGUT), usually of order 10^^ GeV. If <mgut> is set to be msusy, 
the SUSY breaking parameters will be set at Msusy- 

For users that are not familiar with C++, we note that the executable interface allows 
the calculation at just one parameter point in SUSY breaking space. If scans are required. 
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the user can either call SOFTSUSY from a shell script or use a system call from a main 
C program to the executable. Alternatively, a main program showing an example of a scan 
is provided. C++ beginners should note in the following that "method" means function, 
that objects contain a list of data structures and functions and that for a user to access 
(change or reference) the data encoded in an object, one of its functions should be called. 
Such functions are given in tables of the following appendices. 

A.l Input file 

If, as recommended, the SLHA option is used for input, the user may add a SDFTSUSY- 
specific block to the input file in the following format, with bracketed entries replaced by 
double precision values: 

Block SOFTSUSY # SOFTSUSY specific inputs 

1 <TOLERANCE> # desired fractional accuracy in output 

2 <MIXING> # quark mixing option 

3 <PRINTOUT> # gives additional verbose output during calculation 

4 <QEWSB> # change electroweak symmetry breaking scale 

5 <INCLUDE_2_L00P_SCALAR_C0RRECTIDNS> # Full 2-loop running in RGEs 

The fractional numerical precision on masses and couplings output by SOFTSUSY is bet- 
ter than TOLERANCE, which sets the accuracy of the whole calculation. The iteration 
of each physical SUSY particle mass is required to converge to a fractional accuracy 
smaller than TOLERANCE. Sub-iterations are required to converge to a better accuracy than 
10~^ X TOLERANCE for sw and 10~^ x TOLERANCE for fi. The accuracy of the Runge-Kutta RGE 
changes from iteration to iteration but is proportional to the value of TOLERANCE. Values 
between 10^^ and 10^^ are common, lower values mean that SOFTSUSY takes significantly 
longer to perform the calculation. 

The next parameter MIXING determines what Mz boundary condition will be used for 
the quark Yukawa matrix parameters. MIXING=-1.0 sets all Yukawa couplings to zero at Mz 
except for the third-family ones (dominant third-family approximation). MIXING=0.0 sets 
the quark mixings to zero but includes the first two family's diagonal terms. MIXING=1. 0,2.0 
sets all the mixing at AIz to be in the up-quark or down-quark sector respectively, as in 
eq. (pD. 

Setting PRINTOUT to a non-zero value gives additional information on each successive 
iteration. If PRINTOUT>0, a warning flag is produced when the overall iteration finishes. 
The predicted values of M^°'^ and tan /3{Ms) after iteration convergence are also output^. 
The level of convergence, fi{Ms), m'^{Ms) and Mz are output with each iteration, as well 
as a flag if the object becomes non-perturbative. PRINT0UT>1 produces output on the fine- 
tuning calculation. The predicted values of M^°'^ and tan j3{Ms) are output with each 
variation in the initial inputs. A warning flag is produced when a negative-mass squared 
scalar is present. PRINT0UT>2 prints output on the sub-iterations that determine fi{Ms) 
and swiMs), and flags the nature of any tachyons encountered. Values PRINTOUT>0 are 
only required if additional diagnostics are required for debugging purposes. 



^Note that the input value of tan /3 is the value at Mz- 
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QEWSB may be used to multiplicatively change the scale Mg at which the Higgs potential 
is minimised and sparticle masses calculated, or it may alternatively be used as a fixed 
scale. This can be useful if one wants to examine the scale dependence of the results pO| . 
Setting qEWSB=3; < Mz/l GeV, sets Ms = x^m^^ {Ms)m^^ (Ms). Values from 0.5 to 2 are 
common. If Ms < Mz results from the above expression, Ms = Mz is used. If on the 
other hand a fixed scale is desired from Ms, setting QEWSB > Mz/l GeV in the input results 
in Ms =qEWSB GeV being fixed. 

If INCLUDE_2_L00P_SCALAR_CQRRECTI0NS is switched off (o) as in the default case, 2-loop 
RGEs are used for the Higgs and gaugino masses, n, Yukawa and gauge couplings but 
1-loop RGEs are used for other MSSM parameters. Switching on the 2-loop corrections (l) 
results in a full 2-loop RGE evolution, but slows the calculation by a factor of approximately 
three. 

If softpoint.x is used without the SLHA interface, default Standard Model inputs are 
used from the files def .h and lowe.h. The low energy data is encoded in a QedQcd object 
and must be provided. The default numbers supplied and contained in the QedQcd object 
are are given in units of GeV and running masses are in the MS scheme. For the bottom 
and top masses, either the running mass or the pole mass must be supplied as an input. 
The type of mass not given for input is calculated by SOFTSUSY at the 3-loop QCD level. 



We recommend, along with ref. |32|, that the running mass be used for ruf, since there are 
smaller theoretical errors in the extraction of this quantity from experiment than the pole 
mass. The scale dependent quantities in this object are then evolved to Mz by the method 
toMz, to provide the low-scale empirical boundary condition for the rest of the calculation. 
The MSSM spectrum calculated depends most crucially upon Mz,a{Mz),ctsiMz) and 
input third family fermion masses. 



B. Sample Output 

For the recommended SLHA option, the conventions for the output are fully explained 
in Ref. [^]. We therefore present the non-SLHA compliant SOFTSUSY output for a single 
example calculation which can be run by the command 

./softpoint.x sugra 100 250 -100 10 unified 1 



The parameters here have been chosen from SUGRA point SPSla |33|. The output ob- 
tained was 



SOFTSUSY SUGRA calculation 
mgut = 2.393247e+16 



Gravitino mass M3/2: O.OOOOOOe+00 
Msusy: 4.667356e+02 MW: 8.039821e+01 
Data set: 



mU: 1.390102e-03 

mD: 2.764078e-03 

mE: 5.026670e-04 

aE: 1.279250e+02 



mC: 6.353475e-01 mt 

mS: 6.051876e-02 mB 

mM: 1.039356e-01 mT 

aS: 1.176000e-01 Q 



1.638989e+02 mt'pole: 1.714000e+02 

2.884361e+00 mb(mb): 4.200000e+00 

1.751633e+00 mb'pole: 4.938977e+00 

9.118760e+01 mT'pole: 1.776990e+00 
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loops: 3 thresholds: 1 



SUSY breaking MSSM parameters at Q: 9.118760e+01 

UA(3,3) : 

O.OOOOOOe+00 O.OOOOOOe+00 
2.711465e+00 O.OOOOOOe+00 
O.OOOOOOe+00 -4.956895e+02 



-5.932554e-03 
O.OOOOOOe+00 
O.OOOOOOe+00 
UD(3,3) : 

-1.491378e-01 
O.OOOOOOe+00 
O.OOOOOOe+00 
UE(3,3) : 

-7.385460e-03 
O.OOOOOOe+00 
O.OOOOOOe+00 
iiiQLsq(3,3) : 
3.494430e+05 
O.OOOOOOe+00 
O.OOOOOOe+00 
mURsq(3,3) : 
3.298251e+05 
O.OOOOOOe+00 
O.OOOOOOe+00 
mDRsq(3,3) : 
3 . 276634e+05 
O.OOOOOOe+00 
O.OOOOOOe+00 
mLLsq(3,3) : 
3.902935e+04 
O.OOOOOOe+OO 
O.OOOOOOe+00 
mSEsq(3,3) : 
1.857551e+04 
O.OOOOOOe+00 
O.OOOOOOe+00 



O.OOOOOOe+00 
-3.265318e+00 
O.OOOOOOe+00 

O.OOOOOOe+00 
-1.527041e+00 
O.OOOOOOe+00 

O.OOOOOOe+00 
3.494406e+05 
O.OOOOOOe+00 

O.OOOOOOe+00 
3.298228e+05 
O.OOOOOOe+00 

O.OOOOOOe+00 

3.2766086+05 
O.OOOOOOe+00 

O.OOOOOOe+00 

3.902814e+04 
O.OOOOOOe+00 



O.OOOOOOe+00 
O.OOOOOOe+00 
-1.296751e+02 

O.OOOOOOe+00 
O.OOOOOOe+00 
-2.624597e+01 

O.OOOOOOe+00 
O.OOOOOOe+00 
2.900203e+05 

O.OOOOOOe+00 
O.OOOOOOe+00 
2.126351e+05 

O.OOOOOOe+00 

O.OOOOOOe+00 
3.232042e+05 

O.OOOOOOe+00 

O.OOOOOOe+00 
3.866974e+04 



O.OOOOOOe+00 O.OOOOOOe+00 

1.857307e+04 O.OOOOOOe+00 
O.OOOOOOe+00 1.784658e+04 
m3sq: 2.434581e+04 mHlsq: 3 . 148341e+04 inH2sq: 
Gaugino masses(l,3): 

9.815767e+01 1.885414e+02 6.307201e+02 
Supersymmetric paramieters at Q: 9.118760e+01 
Y"U(3,3) : 

O.OOOOOOe+00 O.OOOOOOe+00 
3.620574e-03 O.OOOOOOe+00 
O.OOOOOOe+00 9.027082e-01 



-1.555013e+05 



7.9215956-06 
O.OOOOOOe+00 
O.OOOOOOe+00 
Y"D(3,3) : 
1.5751306-04 
O.OOOOOOe+00 
O.OOOOOOe+00 
Y"E(3,3) : 



0.0000006+00 
3.4487056-03 
O.OOOOOOe+00 



O.OOOOOOe+OO 
O.OOOOOOe+00 
1.479975e-01 
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2.864485e-05 
O.OOOOOOe+00 
O.OOOOOOe+00 



O.OOOOOOe+00 
5.922845e-03 
O.OOOOOOe+00 



O.OOOOOOe+00 
O.OOOOOOe+00 
1.025530e-01 



higgs VEV: 2.494075e+02 tan beta: l.OOOOOOe+01 smu: 3.464099e+02 
gl: 4.592134e-01 g2 : 6.430272e-01 g3: 1 . 139490e+00 
thresholds: 3 #loops: 2 

Physical MSSM parameters : 

mh'O: 1.095488e+02 mA~0: 3.926544e+02 mH~0: 3.936189e+02 mH~+-: 4.011008e+02 
alpha: -1 . 142202e-01 
sneutrinos (1 , 3) : 

1.861463e+02 1.861434e+02 
mU~(2,3) : 
5.634127e+02 
5.465897e+02 
mD~(2,3) : 
5.688527e+02 
5.464146e+02 
mE~(2,3) : 
2.023270e+02 
1.440522e+02 
thetat : 
mGluino 



5.634110e+02 
5.465880e+02 

5.688511e+02 
5.464128e+02 



1.852040e+02 

5.858263e+02 
3.999008e+02 

5.140942e+02 
5.462770e+02 



2.023417e+02 
1.440445e+02 



2.062159e+02 
1.345834e+02 
5.929027e-01 thetab: 3.216005e-01 thetatau: 
6.066386e+02 



-2.855378e-01 



charginos (1 , 2) : 

1.795193e+02 3.762019e+02 
thetaL: -4 . 203557e-01 thetaR: -2 . 429011e-01 
neutralinos (1 ,4) : 

9.691276e+01 1.797635e+02 -3 . 562353e+02 3.743742e+02 
neutralino mixing matrix (4,4): 

9.848209e-01 1.091056e-01 -6 . 158751e-02 -1 . 201277e-01 
-5.833843e-02 9.386643e-01 9.171542e-02 3.272526e-01 

1.532705e-01 -2 . 827469e-01 6.943969e-01 6.437199e-01 
-5.685500e-02 1.644979e-01 7.110613e-01 -6 . 812487e-01 
Isp is neutralino of mass 9.691276e+01 GeV 



Firstly, the output details the input parameters, starting with MIXING and TOLERANCE. 
After the output of the input QedQcd object and the value of Mqut, the result of the 
iteration algorithm in sec. ^is output in the form of a MssmSoftsusy object, unless the user 
is using the SLHA option (in which case, the output will be in SLHA format). The soft 
SUSY breaking parameters were defined in sec. 2.2, and are listed in appendix p.7| . First of 
all, the soft SUSY breaking parameters are displayed. In order, they are the up, down and 
charged lepton trilinear scalar matrices (in units of GeV). Next come the mass squared 
values of the left-handed squarks, right-handed up squarks, right-handed down squarks, 
left-handed sleptons, right-handed charged sleptons in GeV^. m|, m|^^, m^^^ and gaugino 
mass parameters follow. The parameter, 7713/2 (not used here) is the VEV of a compensator 
superfield in anomaly-mediation and completes the SUSY breaking parameter list. 
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Supersymmetric parameters (see sections 2J,E^) are displayed next: Yukawa matri- 
ces , tan/J,(7j, the accuracy level of the calculation, bilinear superpotential fj, 
parameter, renormalisation scale and maximum number of loops used for RGE. 

Physical MSSM parameters follow. The pole masses and mixing parameters are pre- 
viously listed in sec. and are detailed in appendix [E.7| . All masses are in units of GeV, 
and all mixing angles are given in radians. Respectively, there is: rrifiO, ruj^o, rrifjo, rrijjo 
and a. Scalar sparticle masses m^;, mn, mj, mg follow, as well as the mixing angles 9t, Oh, 
9r- The gauginos are listed (in order): m^, m^±, 9l, Or, m^o and O. The DR Higgs VEV 
v{Ms) is then listed, followed by the MS low energy data used as a boundary condition at 
Mz- Finally, the identity of the lightest supersymmetric particle is shown, together with 
its mass. 

Problems (if any) with the parameter point are printed in the final line (none in this 
case, indicating a viable point). 

B.l Problem flags 

Any associated problems such as negative mass-squared scalars or inconsistent EWSB are 
flagged at the end of the output. None of these are printed because the SPSla point 
displayed has none of these problems. We now list the problems, indicating their meaning: 

• If No convergence appears, then SOFTSUSY is indicating that it didn't achieve the 
accuracy of TOLERANCE within less than 40 iterations. The output of the code is 
therefore to be considered unreliable and it is not clear from the output whether 
the point is allowed or disallowed, despite the presence or absence of other warning 
messages. This error flag often appears near the boundary of electroweak symmetry 
breaking, (where h[Msusy) = 0)), where the iterative algorithm is not stable. To 
calculate the position of the electroweak symmetry boundary, one should interpolate 
between regions a small distance away from it. 

• Non-perturbative indicates that SOFTSUSY encountered couplings reaching Landau 
poles when evolving, and could not calculate any further. Any results obtained using 
perturbation theory (for example those of SOFTSUSY) therefore cannot be trusted. 

• Infra-red quasi fixed point breached indicates that the parameter point is at a 
Landau pole of a Yukawa coupling. This should not be a problem provided no other 
errors are flagged. 

• muSqWrongSign indicates that the Higgs minimisation conditions imply that < 0, 
meaning that the desired electroweak minimum is not present in the model. The 
model is ruled out. 

• The t achy on flag occurs when a scalar particle other than the Higgs has acquired a 
negative mass squared, when < or when a pole Higgs masses is imaginary. The 
model is ruled out. 

• noRhoConvergence is flagged when SOFTSUSY cannot calculate the p parameter and 
determine gauge couplings from data, typically because of tachyons or infinities that 
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have crept into the calculation. The other problems are serious enough to rule the 
model out. 



higgsUf b and bProblem indicate that the desired electroweak minimum is in fact a 
saddle point of the potential, thus the model is ruled out. 

mgutOutOf Bounds is flagged if the value of the gauge unification scale predicted by 
eq. |1| is outside the range 10'' GeV < Mx < 5 x 10^'^ GeV. The GUT-scale has 



been set to the appropriate limit, and the SOFTSUSY numbers cannot be trusted. 

Thus flags other than No convergence. Infra-red quasi fixed point breached or Non-perturbative 

indicate an unphysical minimum of the scalar potential, effectively ruling the model point 

out. 



C. Sample Program 

We now present the sample program from which it is possible to run SOFTSUSY in a simple 
fashion. The program we presents here performs a scan in the variable tan with other 
parameters as in mSUGRA point SPSla. It then prints the four pole Higgs masses as a 
function of tan /3 in the standard output channel. If there are any problems with the param- 
eter point, the program prints out these instead of the Higgs masses. The most important 
features of the objects are described in appendix ^ The sample program main . cpp has the 
following form: 

#include <iostream> 
#include "mycomplex.h" 
#include "def.h" 
#include "linalg.h" 
#include "lowe.h" 
#include "rge.h" 
#include "softsusy.h" 
#include "softpars.h" 
#include "susy.h" 
#include "utils.h" 
#include "numerics .h" 



/// global variable declaration 
namespace softsusy { 

/// no quark mixing and no verbose output 

int MIXING = 0, PRINTOUT = 0; 

/// fractional accuracy required 

double TOLERANCE = l.Oe-3; 

/// decay constant of muon 

double GMU = 1.16637e-5; 

/// there are two possible conventions: if QEWSB > MZ, its value is assumed 
/// in GeV and used as a constant MSUSY. Otherwise, it MULTIPLIES the usual 
/// MSUSY value, of root(mstopl mstop2) 
double QEWSB = 1.0; 
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/// Do we include 2-loop RGEs of *all* scalar masses and A-terms, or only the 
/// scalar mass Higgs parameters? (Other quantities all 2-loop anyway) : the 
/// default in SOFTSUSY 3 is to include all 2-loop terms 
bool INCLUDE_2_L00P_SCALAR_C0RRECTI0NS = true; 

/// number of loops used to calculate Higgs mass and tadpoles. They should be 

/// identical for a consistent calculation 

int numHiggsMassLoops = 2, numRewsbLoops = 2; 

/// global pole mass of MZ in GeV - MZCENT is defined in def .h 

double MZ = MZCENT; 

/// end of global variable declaration 



int mainO -[ 
/// Sets format of output: 6 decimal places 
outputCharacteristics(6) ; 

cerr « "SOFTSUSY" « SOFTSUSY_ VERSION 

« " test program, Ben Allanach 2002\n"; 
cerr « "If you use SOFTSUSY, please refer to B.C. Allanach, \n" ; 
cerr « "Comput. Phys . Commun. 143 (2002) 305, hep-ph/0104145\n" ; 
cerr « "B.C. Allanach and M. Bernhardt, arxiv: 0903. 1805 [hep-ph] \n\n" ; 

/// Parameters used: mSUGRA parameters 

double ml2 = 250., aO = -100., mGutGuess = 2.0el6, tanb = 10.0, mO = 100.; 
int sgnMu = 1; ///< sign of mu parameter 

int numPoints = 10; ///< number of scan points 

QedQcd oneset; ///< See "lowe.h" for default definitions paraimeters 

/// most important Standard Model inputs: you may change these and recompile 
double alphasMZ = 0.1187, mtop = 173.4, mbmb =4.2; 
oneset. setAlpha(ALPHAS, alphasMZ) ; 
oneset . setPoleMt (mtop) ; 
oneset . setMass (mBottom, mbmb) ; 

oneset .toMzO ; ///< Runs SM fermion masses to MZ 

/// Print out the SM data being used, as well as quark mixing assumption and 
/// the numerical accuracy of the solution 

cout « "# Low energy data in SOFTSUSY: MIXING=" « MIXING « " TOLERANCE=" 
« TOLERANCE « endl « oneset « endl; 

/// Print out header line 

cout « "# tan beta mh mA mHO mH+-\n"; 

int i; 

/// Set limits of tan beta scan 

double startTanb =3.0, endTanb = 50.0; 
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/// Cycle through different points in the scan 
for (i = 0; i<=numPoints ; i++) { 



tanb = (endTanb - startTanb) / double (numPoints) * doubled) + 
startTaaib; // set tan beta ready for the scan. 

/// Preparation for calculation: set up object and input parameters 
MssmSof tsusy r; 
DoubleVector pars (3); 

pars(l) = mO; pars (2) = ml2; pars (3) = aO; 

bool uni = true; // MGUT defined by gl(MGUT)=g2(MGUT) 

/// Calculate the spectrum 

r . lowDrg(sugraBcs , mGutGuess, pars, sgnMu, tanb, oneset, uni); 

/// check the point in question is problem free: if so print the output 
if ( ! r . displayProblemO . test ) 

cout « tanb « " " « r .displayPhysO .mhO « " " 

« r .displayPhysO .mAO « " " 

« r .displayPhysO .mHO « " " 

<< r .displayPhysO .mHpm << endl; 

else 

/// print out what the problem(s) is(are) 

cout « tanb « " " « r . displayProblemO « endl; 

} 

> 



After an initial introductory print-out, the variables specifying the super symmetry breaking 
parameters are specified. For these, the same notation as appendix |A| is used. Next, the 
important Standard Model inputs are defined and combined with the defaults already 
present in the QedQcd object. The top running mass is calculated from the pole mass and 
Standard Model fermion masses and gauge couplings are then run up to Mz with the 
method toMz. 

If gaugeUnif ication=true, sof tsusy will determine mGutGuess 
from electroweak gauge unification, using the mGutGuess value 
supplied as an initial guess. The user can supply a void func- 
tion that sets the supersymmetry breaking parameters from 
an input DoubleVector. In the sample code given above, this 
function is sugraBcs and is applied to the MssmSoftsusy ob- 
ject at the scale Mgut, which will be determined by the 
user (but mGutGuess will be used as an initial guess). Other 
examples of available boundary conditions are given in Ta- 
ble The user must supply a DoubleVector containing the 
numerical values of the arguments, correctly ordered as in 
Table |2[ sugraBcs, for example, calls the MssmSoftsusy method standardSugra(mO , ml2, 
aO), which sets all scalar masses equal to mO, all gaugino masses to ml2 and all trilinear 



name 
sugraBcs 
amsbBcs 
gmsbBcs 



arguments 
mo, mi/2, Ao 

f^mess, ^ 



Table 2: SUSY breaking 
boundary conditions available 
to the user, detailing argu- 
ments in order. The asterisk 
denotes additional information 
in the text. 
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scalar couplings to aO, in the standard universal fashion. If the user desires to write his or 

her own boundary condition, it must conform to the prototype 

void userDef inedBcs (MssmSof tsusy & m, const DoubleVector & inputs) 

The method lowOrg drives the calculation after, which tan/? and the Higgs masses are 

printed out. 

D. Switches and Constants 

The file def .h contains the switches and constants. If they are changed, the code must be 
recompiled in order to use the new values. Table |^ shows the most important parameters in 
def .h, detailing the default values that the constants have. and Mz have been obtained 
using the latest particle data group numbers [p^]. 



variable default description 



ARRAY_BOUNDS_CHECKING off Vector and Matrix bounds checking 

EPSTOL 10^^^ Underflow accuracy 

GMU 1.16637 10^^ G^, Fermi constant from muon decay 

MZCENT 91.1876 Pole mass of the Z° boson Mz- 

Table 3: Switches and constants. Starred entries have more explanation in the text. is in units 
of GeV2 and Mz in GeV. 



E. Object Structure 

We now go on to sketch the objects and their relationship. This is necessary information 
for generalisation beyond the MSSM. Only methods and data which are deemed important 
for prospective users are mentioned here, but there are many others within the code itself. 

E.l Linear Algebra 

The SOFTSUSY program comes with its own linear algebra classes: Complex, DoubleVector, 
DoubleMatrix, ComplexVector, ComplexMatrix. Constructors of the latter four objects in- 
volve the dimensions of the object, which start at 1. Complex objects are constructed with 
their real and imaginary parts respectively. For example, to define a vector aj=i^2,3; a 
matrix m-j=i...3j=i...4 of type double and a Complex number h = 1 — i: 

DoubleVector a(3) ; 
DoubleMatrix m (3, 4); 
Complex b(1.0, -1.0); 

Obvious algebraic operators between these classes (such as multiplication, addition, sub- 
traction) are defined with overloaded operators *, +, - respectively. Elements of the vector 
and matrix classes are referred to with brackets (). DoubleVector and DoubleMatrix classes 
are contained within each of the higher level objects that we now describe. 

E.2 General Structure 
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From a RGE point of view, a particu- 
lar quantum field theory model consists of a 
set of couplings and masses defined at some 
renormalisation scale fi. A set of /3 func- 
tions describes the evolution of the parame- 
ters and masses to a different scale fi' . This 
concept is embodied in an abstract RGE ob- 
ject, which contains the methods required 
to run objects of derived classes to different 
renormalisation scales. The other objects 
displayed in figure ^ are particular instances 
of RGE, and therefore inherit from it. QedQcd 
objects consist of data on the quark and lep- 
ton masses and gauge couplings. They con- 
tain the P functions for running in an effec- 
Figure 2: Heuristic high-level object structure 

of SOFTSUSY. Inheritance is displayed by the ^^^^ QEDxQCD theory below m*. An ob- 
jjjjgg ject of class MssmSusy contains the Yukawa 

couplings, and the three gauge couplings of 
the MSSM. It also contains the superpotential /i term (not to be confused with the renor- 
malisation scale), tan /?, the ratio of the two Higgs doublet VEVs as well as f = \/v'f + v^- 
Its (3 functions are valid in the exact SUSY limit of the MSSM. The major part of the 
code resides within the MssmSof tsusy class. Objects of this type have all the function- 
ality of MssmSusy, with soft SUSY breaking terms and theoretical boundary conditions 
contained in the inherited class Sof tParsMssm. It also contains an object of type QedQcd 
which contains weak scale empirical data. Code in the MssmSof tsusy class organises and 
performs the main part of the calculation. MssmSof tsusyAltEwsb objects are a slight vari- 
ant of MssmSof tsusy: it takes non-universal Higgs mass boundary conditions at the SUSY 
breaking scale. FlavourMssmSof tsusy objects are also slight variants: they have full flavour- 
mixed output and input. In the following, we provide basic information on classes so that 
users may program using SOFTSUSY. Highly detailed and technical documentation on the 
program may be obtained from the SOFTSUSY website. 

E.3 RGE Class 

The data and important methods in RGE are presented in table ^. Each of the higher level 
objects described in this appendix have explicitly named display and set methods that 
are used to access or change the data contained within each object. In table |^ (as in the 
following tables in this section), these accessing methods are listed on the same row as the 
relevant data variable. 

The RGE method runto (mup , eps) will automatically run any derived object to the scale 
mup with a fractional accuracy of evolution eps. In order to define this evolution, any object 
that inherits from an RGE must contain three methods: display, set, beta shown in table |^. 
DoubleVector display () const must return a vector containing all masses and couplings 
of the object, in some arbitrary user-defined order, void set (const DoubleVector & v) 



RGE 



QedQcd 





MssmSusy 












Sof tParsMssm 






MssmSof tsusy 



FlavourMssmSof tsusy 



MssmSof tsusyAltEwsb 
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must set these couplings given a DoubleVector v defined in the same order as the display 
function. DoubleVector betaO const must then return the /5 functions in a DoubleVector 
defined as 

ft = (E.l) 

where denotes any mass or coupling of the model. The ordering of the must be 
identical in each of the three methods. 



data variable 




methods 


double mu= fj, 


renormalisation scale 


setMu 




(GeV) 


displayMu 


int numpars 


number of scale dependent 


setPars 




parameters 


howMany 


int loops 


accuracy of RGE 


setLoops 






displayLoops 


int thresholds 


accuracy level of threshold 


setThresholds 




computation 


display Thresholds 


method 


function 




DoubleVector display () 


displays all running parameters (*) 


void set (DoubleVector) 


sets all running parameters (*) 


DoubleVector beta 


displays beta functions of all 


running parameters (*) 


runto 


runs object to new value of mu 



Table 4: Abstract RGE class. (*) indicates that derived objects must contain these methods 
(see text). 



E.4 QedQcd Class 

The QedQcd class contains a DoubleVector of quark and lepton MS'masses (my = mu^d,e,c,s,ii,t,b,T{fJ')), 
as shown in table ^. Its contents may be printed to standard output or read from standard 
input (with the same format in each case) by using the operators « or », as can all the 
non-abstract objects mentioned in this section. The methods toMzO, toMtO act on an 
initial object defined with each fermion mass m/ defined at a scale 

Q' = max(l GeV,m/(m/)) (E.2) 

and gauge couplings at Mz- 
E.5 MssmSusy Class 

The operators «, >> have been overloaded to write or read a MssmSusy object to/from a file 
stream. Table ^ shows the data variables and important methods contained in the class. 
For the Yukawa and gauge couplings, methods exist to either set (or display) one element 
or a whole matrix or vector of them. 
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data variable methods 



UouDlevector a 


MS gauge couplings 


setAlpha 






displayAlpha 


DoubleVector m 


running fermion masses 


setMass 


m/(/i) 


vector (1. . . 9) (GeV) 


displayMass 


double mtPole, mbPole 


pole top/bottom/tau 


setPoleMt, setPoleMb 


double mtauPole 


mass 


setPoleMtau 


pole pole pole 

m J , , rrVr 


(GeV) 


displayPoleMt displayPoleMb 



displayPoleMtau 



method 



function 



runGauge 
toMt, toMZ 



runs gauge couplings only 

runs fermion masses and gauge couplings 

from Q' to m^°^'^ or Mz 



Table 5: QedQcd class. Q' is defined in the text. 



E.6 SoftParsMssm Class 

The operators «, » have been overloaded to write or read a sof tParsMssm object to/from 
a file stream. Table |^ shows the data variables and important methods contained in the 
class. addAmsbO adds anomaly mediated supersymmetry breaking terms to the model's 
soft parameters. Such terms are proportional to the VEV of a compensator superfield, so 
771,3/2 in table must have been set before addAmsb is used. minimalGmsb(int n5, double 
lambda, double mMess) applies the messenger scale mMess boundary conditions to the soft 
masses in minimal gauge-mediated supersymmetry breaking ||3^. n5 denotes the number 



of 5 © 5 messenger fields that are present and lainbda(A) is as described in ref. [38| 



E.7 MssmSoftsusy Class 

MssmSof tSusy objects contain a structure sPhysical encapsulating the physical information 
on the superparticles (pole masses and physical mixings), as shown in table |8[ Another 
structure of type drBarPars inherits from sPhysical but instead contains information on 
DR masses and mixing angles. MssmSoftSusy objects also contain one of these structures 
for calculational convenience: the information is used in order to calculate loop corrections 
to various masses. As table |^ shows, a method mpzCharginos returns the 2 by 2 com- 
plex diagonalisation matrices U, V that result in positive DR chargino masses, as defined in 
ref. [^. The method mpzNeutralinos is present in order to convert O to the complex matrix 
N defined in ref. that would produce only positive DR neutralino masses. This infor- 
mation, as well as DR third family fermion masses are stored in the drBarPars structure. 
Another structure within MssmSoftsusy of type sProblem flags various potential problems 
with the object, for example the lack of radiative EWSB or negative mass squared scalars 
(excluding the Higgs mass squared parameters). This structure is shown in table [l^. In 
addition, the method test prints out if any of the possible data variables flagging problems 
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data variable 




methods 


DoubleMatrix u, d, e 


Yukawa couplings 


setYukawaElement 




(3 by 3 matrix) 


setYukawaMatrix 






displayVukawaElement 






display YukawaMatrix 


DoubleVector g 


TV r TV /T 1 " 

MbbM gauge couplings 


setAllGauge 


9% 


(1 ... 3) vector 


setGaugeCoupling 






displayGauge 






displayCaugeCoupling 


double smu 


bilinear Higgs superpotential 


setSusyMu 


A* 


parameter 


displaySusyMu 


double tanb 


ratio of Higgs VEVs (at 


setTanb 


tan /3 


current renormalisation scale) 


display! anb 


double hVev 


Higgs VEV 


setHvev 


V 




displayHvev 



method function 

setDiagYukawas calculates and sets all diagonal Yukawa couplings 

given fermion masses and a Higgs VEV 

getMasses calculates quark and lepton masses from Yukawa 

couplings 

getQuarkMixing mixes quark Yukawa couplings from mass to weak basis 

getQuarkMixedYukawas sets all entries of quark Yukawa couplings given fermion 

masses, Higgs VEV and CKM matrix 

Table 6: MssmSusy class. 



are true. The higgsUf b flag is true if 

m\^^2l? ^m\^-2\ml\<^ (E.3) 

is not satisfied, implying that the desired electroweak minimum is either a maximum or a 
saddle-point of the tree-level Higgs potential IQ]. The contents of sPhysical and sProblem 
can be output with overloaded « operators. noConvergence means that the desired accu- 
racy was not reached. nonPerturbative or irqfp flags the existence of a Landau pole in the 
renormalisation group evolution, and the calculation is not perturbatively reliable so any 
results should be discarded. All other problems except noConvergence and nonPerturbative 
should be considered as grounds for ruling the model out. noRhoConvergence occurs when 
the pseudo-scalar Higgs has a negative mass squared (i.e. an invalid electroweak vac- 
uum). 

MssmSof tsusy data variables and accessors can be viewed in table ^ and the most 
important high-level methods are displayed in table The operators «, » have been 
overloaded to write or read MssmSoftusy objects or sPhysical structures to/from a file 
stream. The driver routine for the RGE evolution and unification calculation is 
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data 




methods 


double m32 


compensator VEV* 


setM32 


^■3/2 


(GeV) 


displayGravitino 


DoubleVector mGaugino 


(1 ... 3) vector of gaugino 


setGauginoMass 


Mi,2,3 


mass parameters 


displayGaugino 


DoubleMatrix ua,da,ea 


(o by 6) matrix oi triimear 


setTrilinearElement 


Ua^Da^Ea 


soft terms (GeV) 


display TrilinearElement 
disnlavSof t A 


LyVJLI.LJJ.dlCLOX Hi J_i O 


(/^ nv 1 TTiJitrir'PS ot Knit 




mURsq , mDRsq , mLLsq 


SUSY breaking masses 


setSof tMassMatrix 


mSEsq 


(GeV2) 


displaySof tMassSquared 


KJ' ("^D' K«)' 






double m3sq,mHlsq,mH2sq 


Bilinear Higgs parameters 


setMSSquared 




(GeV, GeV2, GeV^) 


setMhlSquared 



setMh2Squared 
displayMSSquared 
displayMhlSquared 
displayMh2Squared 



method function 



standardSugra 


Sets all universal soft terms 


universalScalars 


Sets universal scalar masses 


universalGauginos 


Sets universal gaugino masses 


universalTrilinears 


Sets universal soft breaking trilinear couplings 


addAmsb 


Adds AMSB soft terms to current object* 


minimal Gmsb 


Gauge-mediated soft terms used as boundary conditions* 



Table 7: Sof tParsMssm class data and methods. 



double MssmSof tsusy : : lowOrg 

(void (*boundaryConditioii) (MssmSof tsusy &, const DoubleVector &) , 
double mxGuess, const DoubleVector & pars, int sgnMu, double tanb, const QedQcd & 
oneset, bool gaugeUnif ication) 

The user-supplied boundaryCondition function sets the soft parameters according to the ele- 
ments of the supplied DoubleVector at mxGuess, as discussed in appendixj^. If gaugeUnif icationis 
true, the scale that unifies the electroweak gauge couplings is used and returned by the 
function. If gaugeUnif ication is false, the function simply returns mxGuess. pars contains 
a DoubleVector of soft SUSY breaking parameters to be applied as the theoretical bound- 
ary condition. sgnMu is the sign of the superpotential jj, parameter, tanb is the value of 
tw. f5{Mz) required and oneset contains the Mz scale low energy data. 

The fine tuning (as defined in sec. ^) can be calculated with the method 

DoubleVector MssmSof tsusy :: fineTune (void (*boundaryCondition) 
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data variable 



description 



/i^, ^0, if masses 
vector of m/j-^-^ 3 masses 
vectors of m^± ^, i^x°i-i 4 ^respectively 
gluino mass rrig 

4 by 4 orthogonal neutralino mixing matrix O 
(^L.R chargino mixing angles 
9t^b sparticle mixing angles 
double thetatau, thetaH 9^,0 sparticle and Higgs mixing angles 
DoubleMatrix mu, md, me (2 by 3) matrices of up squark, down squark and 

charged slepton masses 
double tlOVlMs, t20V2Ms tadpoles h/vi and t2/v2 evaluated at Ms 

Table 8: sPliysical structure. Masses are pole masses, and stored in units of GeV. Mixing 
angles are in radian units. 



double mhO , mAO , mHO , mHpm 
DoubleVector msnu 
DoubleVector mch.mneut 
double mGluino 

DoubleMatrix mixNeut 
double thetaL, thetaR 
double thetat, thetab 



data variable 



description 



double mt, mb, mtau 

mt{Q),mb{Q),mr{Q) 




Third family fermion masses 


DoubleVector mnBpmz, 

m^o , m + 


mchBpmz 


Absolute neutralino and chargino masses 
(1. . .4, 1. . . 2) vectors 


ComplexMatrix nBpmz 




Neutralino mixing matrix 
(4 by 4 complex matrix) 


ComplexMatrix uBpmz, 

u,v 


vBpmz 


Chargino mixing matrices 
(2 by 2 complex matrices) 


name 




function 


mpzNeutralino 
mpzChargino 




Gives mixing matrices required to make 

neutralino masses positive* 

Gives mixing matrices required to make 



chargino masses positive* 



Table 9: drBarPars structure. Masses are in the DR scheme, and stored in units of GeV. 
Mixing angles are in radian units. Functions marked with an asterisk are mentioned in the 
text. 



(MssmSof tsusy &, const DoubleVector &) , const DoubleVector 
& bcPars, double mx) const 

This function should only be applied to an MssmSof tsusy object which has been processed 
by lowOrg. mx is the unification scale and boundaryCondition is the function that sets 
the unification scale soft parameters, as discussed above. In derived objects, the virtual 
method methodBoundaryCondition may be used to set data additional to MssmSof tsusy from 
the boundaryCondition function. The method outputs the fine-tuning of a parameter ai=i...„ 



-28- 



data variable 


flags 


irqfp 


in a region with a Landau pole 


noConvergence 


the main iteration routine doesn't converge 


noRhoConvergence 


the p iterative routine doesn't converge 


tachyon 


a non-Higgs scalar has negative mass squared 


muSqWrongSign 


p? from eq. (|3.11|) negative 


b 


B from eq. (|3.12|) has incorrect sign 


higgsUf b 


eq. (|E.3|) is not satisfied 


nonperturbat i ve 


a Landau pole was reached below the unification scale 


Table 10: sProblem structure. All data variables are boolean values. 


data 


methods 


double mwPred 


pole Mw prediction setMw 


Mw 


(GeV) displayMw 


double msusy 


Minimisation scale displayMsusy 


Ms 


(GeV) setMsusy 


QedQcd dataset 


Mz boundary condition on setData 
Standard Model couplings displayDataSet 



sProblem problem problem flags displayProblem 



flaglrqfp, f lagB 

f lagNonperturbative 

f lagTachyon, flagHiggsufb 

f lagNoConvergence 

f lagNoMuConvergence 

f lagNoRhoConvergence 

f lagMusqwrongsign 

f lagAllProblems 

DrBarPars f orLoops DR masses and mixings displayDrBarPars 

setDrBarPars 

sPhysical physpars pole masses and mixings displayPhys 

setPhys 

Table 11: MssmSoftsusy class data and accessor methods. 

in the bcPars(n+3) DoubleVector, with the (n + 1, n + 2, n + element of bcPars being the 
fine-tuning with respect to the Higgs potential parameters (/U and B) and the top Yukawa 
coupling {htj respectively. fineTune is an optional feature. sinSqThetaEf f () returns a 
double number corresponding to a full one- loop calculation of the quantity sin^^^^j. It 
does not contain any 2-loop corrections, and may not be accurate enough for precision 
electroweak fits. 

E.8 MssmSoftsusyAltEwsb Class 

The MssmSoftsusyAltEwsb class, which inherits directly from MssmSoftsusy, adds two private 
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name 



function 



lowOrg 

methodBoundaryCondition 
itLowsof t 

sparticleThresholdCorrections 

physical 
calcDrBarPars 
rewsb 
f ineTune 

getVev 

calcSinthdrbar 

calcMs 
printShort 

printLong 

outputFcncs 

lesHouchesAccordOutput 
SinSqThetaEf f 



Driver routine for whole calculation* 
Boundary condition for derived objects* 
Performs the iteration between Mz and 
unification scale 



DR radiative corrections to Standard Model 
couplings at Mz 

Calculates sparticle pole masses and mixings 
Calculates DR pole masses and mixings 
Sets n, B from EWSB conditions 
Calculates fine-tuning for soft parameters* 
and ht 

Calculates VEV at current scale from Z 
self-energy and gauge couplings 
Calculates at current scale from 
gauge couplings 
Calculates Ms 

short list of important parameters printed out 

to standard output in columns 

long list of important parameters printed out 

to standard output in columns 

prints a list of flavour-changing 5 parameters* 



prints output in SLHA |29| format 



calculates then returns sin^ 6 



Table 12: MssmSof tsusy methods and related functions. Functions marked with an asterisk 
are mentioned in the text. 



data variables: muCond and mAcond, both massive parameters in units of GeV. They fix the 
boundary conditions on fJ-{MsusY) and the pole pseudo-scalar Higgs mass m^(pole). In 
this case, mA{MsusY) is extracted from the input m^(pole) via rn\{MsusY) = m^(pole) + 
IlAAiMsusY)^ where IlAAiMsusY) is the MSSM pseudo-scalar self-energy correction. The 
Higgs mass squared parameters mj^^ ^ are not set at Mqut'- they are instead set at 



MsusY by solving the simultaneous equations |3.11| , |3.12| and using the relationship be- 
tween ml{MsusY) and mA{MsusY)- 



= sin^ f5{m\ + M|) + — (1 - sin^ (5) - sin^ /^cos^ (5— - + -Ml) 

Vi V2 2 

mi = cos2/3(m^ + M|) + — (l-cos^/?)-sin2/?cos2/?— -(^2^iM|), (E.4) 



where all quantities in eq. E.4 are running parameters evaluated at MsusY- This option 



is covered under the SLHA input parameters EXTPAR 23 , 26 p9| ] 



-30- 



data 



methods 



f lavourPhysical fv 



sfermion mixing/mass data 



setFlavourPhys 
displayFlavourPhys 



double thetal2, theta23 
double thetalS, deltaCkm 

012, 023,0is,S 



CKM matrix parameters 



setThetal2 
setThetalS 
setTheta23 
setDelta 
display Thetal2 
display Theta23 
display ThetalS 
displayDelta 



double 


thetaB12, tnetaB23 


PMNS matrix parameters 


setTnetaB12 


double 


thetaB13 




setTnetaB13 




, (723 




setThetaB23 








displayThetaB12 








display ThetaB23 








display ThetaB13 


double 


mNuE, mNuMu, mNuTau 


light neutrino masses 


setMnuE 




I'm ' "^i't 




setMnuNu 








setMnuTau 








displayMnuE 








displayMnuMu 








displayMnuTau 


double 


mcMc 


running charm quark mass 


setMcMc 


mdrnc 


) 




displayMcMc 


double 


md2GeV, mu2GeV, ms2GeV 


quark masses at 2 GeV 


setMd2GeV, setMu2GeV 


md{2), 


m„(2),m<j(2) 




setMs2GeV 








displayMd2GeV 








displayMu2GeV 








displayMs2GeV 


double 


mePole, mmuPole 


pole lepton masses 


setPoleMe, setPoleMmu 



pole pole 



displayPoleMe 
displayPoleMmu 



Table 13: FlavourMssmSof tsusy class data and accessor methods. All angles are measured 
in radians and masses are measured in GeV. 



E.9 FlavourMssmSof tsusy Class 

FlavourMssmSof tsusy objects inherit directly from MssmSof tsusy with the addition of the 
variables shown in Table 13. Some of these constitute the angles of the PMNS and CKM 
matrices in the standard parameterisation of eq. 3.7. The lesHouchesAccordOutput method 
has been overloaded to provide flavour violating input and output in accordance with the 
SLHA2 conventions EO]. Note that SOFTSUSY currently does not contain CP-violating 
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complex phases, despite the inclusion of 5. Thus, when CKM angles are input via the 
SLHA2 in the Wolfenstein parameterisation, the magnitude of the 13 entry is fit to sin ^13 



as in eq. 3^ and 6 is set to zero (for now). The relevant method is 



void FlavourMssmSof tsusy : :setAngles 

(double lambda, double aCkm, double rhobar, double etabar) 

Currently, an identical parameterisation is used for the PMNS matrix that describes lepton 
mixing, except 9ij is replaced by Oij. There is currently no provision for a CP-violating 
phase. Effective light neutrino masses are also included in the object. 

An optional feature intended for studies of flavour-changing neutral currents (FCNCs) 
is the method 

void FlavourMssmSof tsusy : : sCkm 

(DoubleMatrix & deltaULL, DoubleMatrix & deltaURR, DoubleMatrix & deltaULR, 
DoubleMatrix & deltaDLL, DoubleMatrix & deltaDRR, DoubleMatrix & deltaDLR) 
const 

which calculates the parameters {^^'l lr RR)ij calculated in the mass-squared-insertion ap- 
proximation (after a rotation to the super CKM basis), as defined in Ref. 

The relationship between the super CKM basis of the mass matrices (including one- 
loop corrections) and the pole-mass basis as described by SLHA2, is contained in the 
structure f lavourPhysical. We list the relevant data in Table |l^. 

data description 

DoubleMatrix dSqMix 6x6 down squark mixing 

DoubleMatrix uSqMix 6x6 up squark mixing 

DoubleVector msD 6 mass-ordered down squark masses 

DoubleVector msU 6 mass-ordered up squark masses 

DoubleMatrix eSqMix 6x6 charged slepton mixing matrix 

DoubleMatrix nuSqMix 3x3 sneutrino mixing matrix 

DoubleVector msE 6 mass-ordered charged slepton masses 

DoubleVector msNu 3 mass-ordered sneutrino masses 

Table 14: f lavourPhysical structure. All mass parameters are pole masses, stored in units 



of GeV. The mixing matrix definitions exactly coincide with those in the SLHA2 [40| and 
describe the transformation between the loop-corrected super CKM basis and the mass 
basis. 
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